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SUMMARY 

The  heat-transfer  and  pressure-drop  characteristics  of  liq,uid  para- 
hydrogen  flowlr^  through  a  heated  vertical  tube  were  investigated  over 
the  following  range  of  conditions:  inlet  pressure,  30  to  70  pounds  per 
square  inch;  average  heat  flux,  1  Btu  per  square  inch  per  second  maximum 
and  temperatui’e  differential,  50°  to  750°  R  between  the  fluid  and  the 
wall.  As  stored  in  the  Dewar,  the  hydrogen  was  approximately  95  percent 
para.  Two  different  test-section  geometries  of  1-foot  length  were  em¬ 
ployed;  one  had  a  0.625-inch  outside  diameter  and  0.065-inch  wall,  and 
the  other  a  0.375-lnGh  outside  diameter  and  0.031-inch  wall. 

Local  pressure  drop  and  heat-transfer  coefficients  for  boiling  hy¬ 
drogen  are  presented.  The  wall  temperature  profile  along  the  tube  was 
radically  different  from  that  observed  in  convective  heat  transfer  with 
single-phase  fluids.  For  the  larger  heat  fluxes,  a  minimum  wall  temper¬ 
ature  occurred  somewhere  near  the  end  of  the  tube,  and  the  maximum  tem¬ 
perature  occurred  near  the  entrance.  A  vapor  binding  or  dry-wall  condi¬ 
tion  was  apparent  in  some  of  the  hydrogen  runs.  An  empirical  correla¬ 
tion  for  predicting  the  two-phase  heat- transfer  coefficient  is  presented 

Pressure-drop  measurements  showed  that  the  frictional  loss  is  small 
compared  with  the  momentum  loss.  In  fact,  a  one-dimensional  momentum 
analysis  can  predict  the  pressure  drop  quite  accui-ately . 


INTRODUCTION 

Liquid  hydrogen  appears  to  be  a  very  attractive  propellant  either 
in  chemical  systems  in  which  it  is  burned  with  an  oxidant,  or  in  mono¬ 
propellant  systems  in  which  it  is  heated  by  some  device  like  a  nuclear 
reactor.  Its  low  atomic  weight,  high  specific  heat,  and  high  heat  of 
combustion  all  contribute  to  high  specific-impulse  values  for  either  of 
these  types  of  propulsion  systems.  The  undesirable  characteristics  of 
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the  prepellwit  include  a  large  specific  tc^Lme  and  the  Ulffieult  1i4riaiing 
procedures  .associated  with  a  low.itehQ)eratu3:’e '  ciyogen . 

Beyond  the  attractive  propulsive  qualities,  hydrogen  appe^s  as  a 
good  coolant  with  an  appreciable  heat-sink  capability.  Since  many  com¬ 
ponents  of  a  propulsion  system  must  be  cooled,  this  is  a  further  .argu¬ 
ment  for  hydrogen  as  a  propellant. 

To  evaluate  liquid  hydrogen  properly  as  a  coolant,  more  experimental 
Information  Is  required  concerning  the  heat-transport  mechanism  for  both 
subcrltical  and  supercritical  conditions.  The  problem  is  particularly 
acute  for  hydrogen,  because  it  is  generally  used  at  conditions  where  den¬ 
sity  changes  rapidly.  For  this  report,  subcrltical  pressures  from  30  to 
70  pounds  per  squau’e  inch  absolute  emu  a  narrow  range  of  bulk  tempera¬ 
tures  (low  subcooling)  are  included  in  the  experimental  conditions.  The 
heat-transfer  characteristics  of  hydrogen  were  measured  in  an  electri¬ 
cally  heated  vertical  tube.  The  tube  was  instrumented  for  surface  tem¬ 
perature  and  pressure-drop  measurements.  The  temperature  difference  be¬ 
tween  the  wall  and  bulk  was  varied  from  approximately  50°  to  750°  R. 

The  maximum  heat  flux  was  1  Btu  per  square  inch  per  second. 

For  the  range  of  investigation,  film  boiling  of  the  hydrogen  and  a 
resulting  two-phase  flow  phenomenon  occurred.  The  heat-transfer  and 
pressure-drop  results  are  discussed  in  terms  of  a  two-phase  flow  model. 

A  correlation  scheme  is  presented  that  is  similar  to  methods  employed 
in  correlating  two-phase  heat-transfer  processes  with  other  fluids. 

One  of  the  variables  employed  in  the  heat-transfer  correlation  is 
the  Martlnelll  two-phase  parameter  Xtt^  which  is  an  index  of  wall  and 
two-phase  shear  resistance  of  the  two-phase  flow.  Originally  presented 
in  reference  1,  this  parameter  was  used  to  predict  two-phase  pressure  drop 
from  an  estimate  of  the  gas-phase  pressure  drop.  Martlnelll 's  friction 
correlation  was  applied  to  a  gaseous  core  and  a  liquid  annulus  flowing 
isothermally  in  a  pipe.  Appendix  D  shows  that  the  same  two-phase  param¬ 
eter  Xtt  evolves  from  an  analysis  similar  to  Martlnelll 's  for  a  liquid 
core  and  a  gaseous  annulus.  Consequently,  the  application  of  this  param¬ 
eter  to  the  hydrogen  two-phase  flow  is  justified. 

Several  researchers  including  Guerrieri  and  Talty  (ref.  2)  and 
Dengler  and  Addoms  (ref.  3)  have  correlated  two-phase  fluid  heat  transfer 
using  the  Martlnelll  parameter.  A  similar  empirical  approach  Is  applied 
to  the  film-boiling  da.ta  of  liquid  hydrogen  reported  herein. 
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APPARASUS 
Plow  System 

The  test  installation  Is  shown  schematically  in  figuje  1;  Liq_uid 
hydrogenj  stored  in  a  large  Dewar,  was  forced  through  the  tlQv  sy&tem  hy 
external  gaseous  hydrogen  pressure.  The  entrance  piping  to  the  test 
section  and  the  test  section  itself  were  vacuum- Jacketed  to  minimize 
heat  leakage.  After  passing  through  the  test  section,  the  hydrogen  en¬ 
tered  a  steam  heat  exchanger  where  it  was  cong^letely  vaporized,  metered, 
and  then  exhausted  Into  the  atmosphere.  All  elements  of  the  stack  dis¬ 
charge  system  were  thoroughly  grounded  to  prevent  buildup  of  static  elec¬ 
tric  sQ.  charge. 

Plow  rate  and  pressure  level  were  controlled  by  remotely  operated 
valves  upstream  .and  downstream  of  the  test  section. 


Test  Section 

A  sketch  of  the  vertical-tube  heat-transfer  test  section  is  shown 
In  figure  2.  Two  test-section  geometries  were  used;  the  first  was  a 
stainless-steel  tube  of  0.625-inch  outside  diameter  and  0.065-inch  wall, 
and  the  second  was  an  Inconel  tube  of  0.375- inch  outside  diameter  and 
0.031-inch  wall.  The  heated  portions  of  both  geometries  were  12  inches 
long,  ^proximately  5  inches  of  unheated  test  section  preceded  the 
heated  portion.  The  smaller  Inconel  tube  was  more  extensively  instru¬ 
mented  for  local  pressure  measurements.  Because  of  the  importance  of 
pressure  data  in  determining  local  heat-transfer  conditions,  the  overall 
data  from  the  smaller  tube  were  more  amenable  to  analysis.  The  test  sec¬ 
tion  was  electrically  heated  by  the  secondary  of  a  65-kilovolt-ampere 
alternating-current  power  supply  (fig.  1).  The  maximum  rated  output  of 
the  secondary  was  13  volts  at  5000  amperes.  The  test  section  was  elec¬ 
trically  insulated  from  the  adjacent  piping  of  the  flow  apparatus  by  thin 
Teflon  gaskets  and  insulating  grommets  in  the  flange  connections. 


Instrumentation 

The  measuring  stations  and  types  of  measurement  used  for  the  0.375- 
Inch  Inconel  tube  are  shown  schematically  in  figure  2.  The  hydrogen  flow 
rate  was  metered  upstream  of  the  test  section  by  a  Venturi  and  downstream 
of  the  heat  exchanger  by  an  orifice.  The  two  flowmeters  served  as  checks 
on  each  other. 

The  test  sections  were  instrumented  with  copper-constantan  thermo¬ 
couples  silver-soldered  to  the  outside  wall.  Additional  thermocouples 


were  placed  on  the  copper  flanges  to  cheek  heat  leakage  from  th6  test 
section.  All  thermocouples  employed  an  atmospheric  holllng  nitrogen 
cold  junction. 

Pressure  taps  were  located  upstream  and  downstream  of  the  test  sec¬ 
tions  and  along  the  heated  portion  of  only  the  Inconel  tube. 


PROCEDURE 

Operating  Conditions 

The  following  ranges  of  operating  conditions  were  investigated: 
flow  rate,  from  0.03  to  0.19  pound-mass  per  second;  inlet  pressure,  from 
30  to  70  pounds  per  square  inch  absolute;  and  electrical  power  input, 
from  0  to  14  Btu  per  second. 


Data  Recording 

During  a  run,  pressure,  flow  rate,  temperature,  and  electrical 
power  were  recorded  on  tape  by  an  automatic  voltage  digitizer  (ref.  4) 
and  were  available  for  immediate  wrlxe-back  on  an  electric  typewriter. 

The  data  tapes  were  processed  by  a  high-speed  computer,  which  converted 
the  data  bits  5nto  the  desired  parameters  for  further  analysis.  Periodic 
checks  of  the  digitized  recording  were  made  on  conventional  self- 
balancing  potentiometers  and  on  a  multichannel  oscillograph. 


Analysis  of  Data 

Symbols  are  listed  in  appendix  A,  and  equations  and  assumptions  em¬ 
ployed  in  reducing  the  experimental  data  appear  in  appendix  B.  The 
reader  should  become  acquainted  with  appendix  B  as  reference  is  made  to 
definitions  and  equations  therein  throughout  the  DISCU'SSICN  OP  RESULTS. 

The  analysis  of  the  heat  losses  peculiar  to  the  apparatus  is  also 
presented  in  appendix  B.  It  is  sufficient  to  say  here  that  the  heat 
leakage  Into  or  out  of  the  test  section  was  negligible.  The  table  in 
appendix  B  is  included  for  the  reader  who  may  be  interested  in  actual 
numerical  values  of  the  heat  loss. 

The  thermodynamic  and  transport  data  of  gaseous  and  liquid  hydrogen 
used  in  the  computation  of  results  are  presented  as  the  last  figures 
herein,  numbered  8(a)  to  (f).  These  data  are  reproduced  herein  because 
they  are  compiled  from  several  sources,  some  published  (refs.  5  to  9) 
and  some  unpublished.  Also,  since  there  are  insufficient  experimental 
property  and  transport  data  at  the  present  time,  the  authors  extrapolated 
these  data  for  their  own  use.  For  exaagjle,  the  only  reliable  viscosity 
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and  speGific-heat  data  for  liqiuid  hydrogen  pertain  to  1-atmosphere  pres¬ 
sure,  hut  the  data  were  applied  to  conditions  at  several  atmospheres. 

As  more  reliable  data  become  available,  the  reader  may  apply  correction 
ratios  to  these  results  to  obtain  results  consistent  with  the  new  data. 


DISCUSSION  OF  RESULTS 

The  experimental  data  and  some  selected  results  for  the  0.375-inch- 
outside-diameter  Inconel  tube  are  tabulated  in  table  I. 


Unll  Ternperature  Profile" 

Figure  3(a)  is  a  family  of  inside- wall  temperature  profiles  for 
various  flow  rates  and  heat  fluxes.  These  I’epresent  typical  tube  temper 
ature  profiles  from  the  0.375-inch  Inconel  tube,  which  are  similar  in 
characteristics  to  those  obtained  with  the  0.625-inch  staiifLess-steel 
tube  (fig.  3(b)).  At  the  lower  heat  fluxes,  the  profile  is  rather  flat; 
but,  as  the  heat  flux  is  Increased,  the  temperature  profile  takes  on  a 
negative  slope.  Opposite  to  single-phase  heat-transfer  results,  the 
tube  wall  is  actually  hotter  at  the  entrance  than  at  the  exit. 

At  first,  it  was  thoitght  that  para  to  ortho  conversion  at  the  tube 
wall  might  be  Influencing  the  tube  temperature  profile.  However,  an 
analysis  of  the  conversion  reaction,  as  presented  in  appendix  C,  led  to 
the  conclusion  that  no  appreciable  endothermic  reaction  could  occur  in 
the  vicinity  of  the  tube  wall,  even  though  the  wall  material  was  cata¬ 
lytic;  the  conversion  rate  is  negligibly  low.  To  verify  this  conclusion 
the  interior  of  the  stainless-steel  test  section  was  brush-plated  with 
gold.  Within  the  experimental  accuracy,  no  differences  in  temperature 
were  observed. 

A  more  reasonable  explanation  for  the  wall  temperature  profile  can 
be  developed  by  considering  a  simple  one-dimensional  model  of  the  flow. 
As  the  hydrogen  enters  the  test  section,  immediately  a  thin  vapor  film 
is  established  adjacent  to  the  wall.  When  saturation  conditions  are 
reached,  a  more  turbulent  exchange  between  the  nonstable  vapor  film  and 
the  liq[uid  core  takes  place.  While  this  exchange  mechanism  probably 
represents  a  nonequilibrium  condition  because  gas  is  generated  at  the 
wall  and  also  in  the  liquid  core  as  pressvire  is  reduced,  equilibrium 
conditions  were  assumed  in  evaluating  quality.  As  the  mixture  moves 
through  the  heated  tube,  more  gas  is  generated,  and  continuity  consid¬ 
erations  demand  that  the  mixture  velocity  increase,  which  decreases  the 
pressure.  A  pressure  decrease,  lowering  saturation  level,  requires  the 
fluid  to  absorb  a  portion  of  the  latent  heat,  thereby  decreasing  density 
and  increasing  velocity.  As  in  any  turbulent  convective  heat-transfer 
process,  increasing  the  velocity  enhances  the  heat-transfer  mechanism; 
consequently,  the  wall  temperature  drops.  The  vapor  generated  before 
saturation  is  assumed  small  and  is  a  good  approximation  for  small 
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subcoollngj  however^  for  large  subcooling^  the  vapor  generated  must  be 
considered  (see  appendix  B). 

There  are  two  in^jortant  exceptions  to  the  condition  where  the  tube 
temperature  diminished  monotonlcally  from  the  inlet  to  the  exit  of  the 
test  section.  For  the  first  exception,  the  liquid  hydrogen  entered  the 
test  section  subcooled  and  consequently  did  not  reach  saturation  condi¬ 
tions  until  some  station  downstream  of  the  tube  entrance.  Thus,  the 
distance  from  the  heater  entrance  to  the  station  where  saturation  condi¬ 
tions  prevail  is  greater  than  zero  (see  definition  of  saturation  length, 
appendix  B) .  It  was  observed  that  the  tube  wall  temperature  increased 
and  pressure  decreased  until  the  saturation-length  station  was  reached. 
Beyond  this  point  the  tube  temperat'ure  diminished  in  the  manner  previ¬ 
ously  described.  It  is  interesting  to  observe  that,  if  the  hydrogen 
entered  the  test  section  at  saturation  conditions  or  with  quality,  the 
tube  temperatures  in  the  inlet  of  the  heated  section  would  be  lower  than 
if  the  hydrogen  were  a  subcooled  liquid.  (See  runs  18-6  and  18-7,  fig. 
3(a);  also  note  w  and  pj^^^  differences,  table  I.)  One  can  conclude 
that  lowering  the  system  pressure  (lowering  the  subcooling)  would  help 
the  cooling  process  at  the  inlet  if  a  subcooled  fluid  were  present. 

(See  runs  20-1,  20-2,  and  20-3,  fig.  3(a),  and  compare  subcooling  through 
use  of  table  I  and  fig.  8(a).) 

The  second  exception  to  the  monotonlcally  decreasing  wall  tempera¬ 
ture  involves  an  abrupt  increase  in  wall  temperature  beyond  the  midpor¬ 
tion  of  the  tube.  The  effect  was  most  noticeable  with  the  0.625- inch 
tube.  Figure  3(b)  shows  this  wall  temperature  profile.  The  figure  also 
shows  how  this  condition  grew  with  increasing  heat  flux.  The  abrupt 
temperature  change  appeared  to  take  place  when  the  estimated  quality 
reached  approximately  50  percent  or  greater.  A  similar  phenomenon  has 
been  observed  with  water  and  other  fluids  in  boiling  processes  (ref.  3), 
and  it  is  called  "vapor  binding"  or  "dry  wall."  So  much  vapor  is  present 
near  the  heating  surface  that  in  the  turbulent  vapor-liquid  mechanism 
practically  no  liquid  particles  get  near  enough  to  the  wall  to  absorb  the 
heat  of  vaporization.  The  heat  being  transferred  appears  principally  as 
sensible  heating  of  the  vapor  phase,  and  the  heat-transfer  coefficient 
more  closely  approaches  the  single-phase  gaseous  value. 


Pressure  Drop 

The  heated  portion  of  the  Inconel  test  section  was  Instrumented 
with  four  pressure  taps.  The  pressure-drop  profiles  obtained  for  several 
flow  rates  and  heat  fluxes  are  shown  in  figure  4.  The  pressure  drop  per 
unit  length  may  be  obtained,  with  slight  modifications,  from  (ref.  10) 

*  (g)  (4 


(1) 


Because  of  the  ajppreoiahle  difference  between  the  eppe^ent  density  of  the 
fluid  at  the  entrance  and  at  the  exit,  the  pressure  drop  could  be  as¬ 
cribed  to  momentum  change,  and  the  viscous  shear  forces  at  the  wall 
could  be  neglected.  Thus,  the  form  can  be  simplified  to  the  one- 
dimensional  momentum  equation: 

a  -JL-  Au  (2) 

The  simplified  form  predicts  pressure-drop  values  that  are  roughly  equiv¬ 
alent  to  the  experimental  values.  The  dashed  lines  calculated  using  the 
one-dimensional  relation  are  superposed  on  the  experimental  data  of  fig¬ 
ure  4  to  Illustrate  the  relative .agreement  of  analysis  and  experiment. 


Heat  Transfer 

Local  values  of  heat-transfer  coefficient  h  are  plotted  as  a 
function  of  the  test-section  length  for  various  mass-flow  rates  and  heat 
Inputs  in  figure  5.  These  are  typical  data  for  both  test  sections  and 
were  selected  to  show  how  the  coefficient  changes  with  heat  flux  as  well 
as  weight  flow.  These  coefficients  for  the  0.375-inch  tube  are  also 
tabulated  in  table  I.  The  heat- transfer-coefficient  results  for  the 
smaller- diameter  test  section  were  Judged  to  be  the  more  reliable  be¬ 
cause  the  local  pressures  were  known  throughout  the  test  section. 

Several  empirical  correlation  schemes  similar  to  those  appearing  in 
references  2  and  3  were  applied  to  the  two-phase  experimental  heat- 
transfer  results.  In  reference  2  the  Martlnelli  parameter  Xtt  and  a 
modified  Reynolds  number  involving  the  mass  fraction  of  vapor  x  were 
used  in  establishing  the  correlation.  Similar  parameters  were  chosen 
for  the  correlation  of  the  hydrogen  data. 

First,  It  has  to  be  established  that  the  Martinelli  flow  parameter 
was  applicable  to  the  two-phase  hydrogen  flow  model.  In  the  original 
development  of  the  parameter  (ref.  l),  the  fluid  model  consisted  of  a 
gaseous  core  and  a  liquid  annulus  adjacent  to  the  wall.  By  a  similar 
derivation  it  was  apparent  that  the  Martinelli  parameter  was  applicable 
to  the  hydrogen  flow  composed  of  a  liquid  core  and  a  gaseous  annulus 
(see  appendix  D) . 

Several  methods  of  computing  a  modified  Reynolds  number  for  the  two- 
phase  flow  were  tried.  The  fi-uid  properties  were  evaluated  at  various 
temperatures  including  wall,  bulk,  and  film  (the  arithmetic  mean  of  hulk 
and  wall)  conditions.  These  Reynolds  numbers,  along  with  corresponding 
Prandtl  numbers,  were  used  in  calculating  an  estimated  Nusselt  number  by 
the  modified  Dlttus-Boelter  equation: 
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Nu^alc^f  =  0.023  Re0-8pr0.4  (3) 

In  selecting  a  Reynolds  numter  that  would  apply  to  a  broad  range  of 
qualities  (mass-percent  vapor)  ^  the  Reynolds  ntanber  should  be  equivalent 
to  a  gaseous  one  when  the  fluid  approaches  the  single-phase  gaseous  con¬ 
dition.  It  was  found  that  this  requirement  could  be  satisfied  by  defin¬ 
ing  Reynolds  number  in  the  following  way: 

^  Pf.m^avD 


where  the  mean  density  of  the  film  is 

Pf,m  ~  1  -  X 

Pf  Pi 

X  is  quality  (greater  than  zero,  as  discussed  later),  is  the  gaseous 
density  evaluated  at  film  conditions,  Pj  is  the  saturated  liquid  den¬ 
sity,  uav  Is  the  average  velocity  of  the  mixture  of  liquid  and  gaseous 
phases,  and  is  the  viscosity  at  gaseous  film  conditions. 

In  addition,  writing  the  density  as  Pf^,n  modifies  the  film  den¬ 
sity  for  the  liquid  droplets  present.  Perhaps  these  droplets,  diffusing 
toward  the  hot  wall,  are  partially  vaporized  and  then  "bounced"  back  into 
the  core.  At  high  qualities,  apparently  very  small  amounts  of  diffusing 
liquid  are  present  in  the  film}  thus,  the  mean  film  density  approaches 
the  gaseous  value. 

Figure  6  shows  a  ratio  of  the  experimental  Nusselt  number  to  the 
predicted  Nusselt  number  (based  on  the  preceding  definition  of  Reynolds 
number)  plotted  as  a  function  of  X^-t.  The  data  points  represent  local 
conditions  at  four  axial  stations  near  the  midportion  of  the  small- 
diameter  Inconel  test  section.  Such  a  selection  of  axial  stations  obvi¬ 
ated  the  end  effects.  The  data  may  be  represented  by  the  following 
equation; 


(4) 

(5) 


NUoqT 2 - f  , 

Kuexp,f  =  0.611  +  1.93  Xtt,f  ^ 

Equation  (6)  is  a  least-squares  fit  to  the  data  from  L  =5.86  to 
9.4  inches  and  is  valid  for  qualities  from  0.05  to  0.9  over  the  range 
of  pressures  covered  by  this  report.  Tiie  equations  may  be  extended, 
with  reservation,  to  qualities  of  0.01  (see  fig.  7). 

It  should  be  noted  that  the  true  case  is  probably  neither  a  sep¬ 
arate  phase  flow  as  assumed  in  the  derivation  of  Xtt  nor  a  homo¬ 
geneous  flow  as  Implied  by  the  introduction  of  x  in  Pf,in.  Hence, 


E-976 


E-976 


9 


the  incorporation  of  hoth  parameters  'based  upon  these  extreme  cases 
gave  correlation  (6). 

As  quality  increases  (X-tt  decreases),  the  ratio  of  Nusselt  numbers 
begins  to  approximate  unity.  Hovever,  a  careftil  look  at  figure  6  shows 
the  Kusselt  number  ratio  to  be  greater  than  unity  at  the  low  Xtt  val¬ 
ues  .  Apparently  the  presence  of  nonequilibrium  quantities  of  colloidal 
particles  alters  the  single-phase  gaseous  heat- transfer  mechanism.  It 
is  expected  that  the  experimental  values  of  Nusselt  number  would  equal 
the  predicted  values  when  the  flow  is  strictly  single-phase  gaseous,  as 
shown  by  reference  11.  Similar  results  were  obtained  for  helium  in  ref¬ 
erence  12. 


It  is  fairly  obvious  that  the  gaseous  film  density  cannot  be  as¬ 
signed  to  the  Reynolds  number  when  the  fluid  is  approaching  the  com¬ 
pletely  liquid  state.  Somewhere  at  high  values  of  X^t^  quality, 

the  Reynolds  number  should  be  redefined  incorporating  liquid  properties. 
Figure  7  includes  the  families  of  local  data  over  the  tube  length.  The 
axial  stations  included  in  the  correlation  begin  about  0.625  inch  from 
the  inlet  (or  where  satiu-ation  occurs)  and  end  1.6  inches  before  the 
tube  exit.  A  least-squares  fit  to  these  data  yields 


NUexp, f 


NUcalc,f _ 

0.706  +  1.60  -  0.123  X|^^f 


This  equation  and  equation  (6)  are  superimposed  on  the  data  of 
figure  7  for  comparison. 


The  feunilies  may  be  represented  by  equation  (6)  with  the  following 
observations: 


(1)  The  predicted  value  of  Nusselt  number  is  usually  lower  than  the 
experimental  at  the  inlet  and  higher  at  the  exit.  This  tendency  in  the 
Nusselt  numbei’  suggests  that  a  length-diameter-ratio  correction  might 
improve  the  correlation.  Additional  tube  geometries  would  be  required 
to  establish  the  effect  of  the  L/D  term  experimentally.  Consequently, 
no  attempt  to  suggest  an  L/d  correction  is  presented  herein- 

(2)  The  corre-lation  applies  for  moderate  subcooling  at  the  pres¬ 
sures  considerably  below  critical  and  to  lower  degrees  of  subcooling  as 
the  pressure  is  increased.  The  sensitivity  of  the  correlation  to  the 
subcooling  is  much  more  pronounced  near  the  critical  pressure. 

As  was  mentioned  in  the  INTRCDUCTION,  X^t  is  a  measure  of  the 
fluid  resistance  of  the  two-phase  flow.  It  has  been  used  to  correlate 
two-phase  fluid-flow  pressure  drop.  In  its  use  as  a  parameter  for  two- 
phase  heat-transfer  correlations,  it  may  be  interpreted  as  a  correction 
to  the  shear  terms  in  the  Dittus-Boelter  equation. 
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A  design  procedure  is  treated  In  some  detail  In  appendix  B*  This 
procedure,  while  nearly  comprehensive,  must  necessarily  he  modified  hy 
the  designer  to  suit  the  particular  application. 


SUMMARY  OP  RESULTS 


1.  The  tui'hnlent  two-phase  para-hydrogen  heat- transfer  data  are 
correlated  using  the  Martinelli  parameter  Xtt,f  ^ind  the  modified 
Dlttus-Boelter  equation  NU(.J^2.c,f‘  data  are  represented  hy 


Nu, 


Nu, 


exo.f  “  r> 


calc,f 


y.  ,  _ 


where 

Nw,f  =  0-023  Re0-8pr0-4 

and 


(6) 


(3) 


Re 


P-  u  D 
f,m  av 


Pf 


(4) 


Equation  (6)  holds  at  least  over  the  range  of  conditions  studied;  inlet 
pressure,  30  to  70  pounds  per  square  inch  absolute;  average  heat  flux,  1 
Btu  per  square  inch  per  second  maximiun;  and  temperature  differential, 

50°  to  750°  R  hetv;een  fluid  and  wall. 


2.  V/all  temperature  profiles  were  generally  high  at  the  inlet  and 
decreased  toward  the  exit;  this  trend  is  just  the  opposite  of  that  usu¬ 
ally  observed  for  convective  heat  transfer  with  single-phase  fluids. 

3.  Pressure-drop  measurements  indicated  momentum  losses  to  he 
greater  than  frictional  loss.  A  one-dimensional  momentum  analysis  can 
be  used  to  predict  pressure  drop. 

4.  Heat  absorption  due  to  para- ortho  conversion  is  small  and  at 
least  within  the  accuracy  of  the  experiment. 

5.  Appendixes  present  Justification  of  assumptions,  methods  of  cal¬ 
culation,  para- ortho  hydrogen  conversion  analysis,  and  a  design 
procedure. 


Lewis  Research  Center 

National  Aeronautics  and  Space  Administration 
Cleveland,  Ohio,  January  11,  1961 
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Nuoalc,f 
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Nusselt  umber  compuied  froffl,  modified  ti'i-ttusxBaglter  y 

tion  using  film  temperature  to  evaluate  properties 
(eq.  (3)) 

experimental  Nusselt  number,  ^exp^/^f 
number  of  voltage  taps 
Prandtl  number,  Cpn/K 
static  pressure,  Ib/sq  in.  or  mm  Hg 
pressure  gradient,  (ib/sq  in. )/in. 
heat  flow,  Btu/sec 
heat  flux,  Btu/(sq  in.) (sec) 

gas  constant,  82.05  (atm)(cc)/(g-mass)(mole)(°K) 

Reynolds  number  (eq.  (i)) 
conversion  rate,  (g-mass) (mole )/(cin^) (sec) 
tube  radius,  in. 
temperature,  OR  or  °K 
temperature  difference,  °R  or  °K 
time,  sec 
velocity,  ft/sec 
volume ,  cc 

mass-flow  rate,  lb-mass/ sec 
quality,  percent  vapor  by  mass 
distance  from  wall 
form  factor,  liquid  area 
form  factor,  gaseous  area 

correction  coefficient  for  hydraulic  diameter 
boundary-layer  thickness,  cm 
eraisslvity 
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l8o1SteBma:L  fsletlon  faxst  ■  """  . .  '  I  “ 

K  thermal  conductivity,  Btu./(ft) (sec) (9R) 

/<U  thermal  conductivity  of  materlalj,  Btu/(ft)(6ee)('^) 

Ty  heat  of  vaporization,  Btu/lb-maas 

(i  viscosity,  lb-mass/ (ft) (sec) 

p  density,  Ib-mass/cu  ft 

a~  mean  film  density  (eq,.  (s)),  lb-mas s/cu  ft 

I  ,m 

0  Stefan-Boltzmann  constant,  0.173X10“8  Btu/(hr)(sq  ft)(9R^) 

T  half-life,  sec 

Xtt  Martlnelll  two-phase  parameter  (both  phases  turbulent) 

(eq.  (D16)) 


Subscripts ; 
av  average 

b  bulk  temperature 

Cu  copper 

can  vacuum  container  enclosing  test  section 

e  exit 

eq  equilibrium 

exp  experimental 

F  flange 

f  film  conditions,  arithmetic  mean  between  wall  and  bulk 

temperature 

g  gas 

i  inside  wall  of  test  section 

in  inlet 


L 

Z 

o 

or 

ortho 

para 

pipe 

R 

sat 

t 

tot 

tp 

V 

V 

W 

w 

3 

4 

5 


station  along  heated  length 
liquid 

outside  wall  of  test  section 
orifice 

ortho-hydrogen 

para-hydrogen 

test  section 

radiation 

saturation 

time 

total 

two  phase 

Venturi 

voltage  tap 

voltage-tap  wire 

wall 

thermocouple  location  l/2  inch  from  flange 
thermocouple  location  at  flange 
instrumentation  Junction  strip 


.  -  », 
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CCMHJTATIQNAL  JROCBDUBE 
Inside-Wall  Temperature 

The  tube  surface  thermocouples  were  mounted  on  the  exterior  surface. 
Values  of  the  Inside-wall  temperature  were  confuted  from  the  following 
expression  (ref.  13),  which  presupposes  uniform  radial  heat  flux  through 
the  tube  wall  (tube  generating  power) ; 

T.  =  T..  -  7.  -S- 

where 


Saturation  Length 

Before  computing  the  local  density  and  quality  of  hydrogen  at  any 
axial  station  in  the  tube,  some  simplifying  assumptions  are  required  to 
establish  where  vapor  accumulation  begins. 

It  is  recognized  that,  for  the  rather  appreciable  differences  be¬ 
tween  the  wall  and  the  bulk  temperatures,  boiling  probably  begins  at  the 
entrance  of  the  test  section.  However,  if  the  local  bulk  ten^serature  in 
the  tube  entrance  is  less  than  the  saturated  temperature,  subcooled  boil¬ 
ing  occijr.?.  For  most  of  the  actual  experimental  runs,  some  degree  of 
subcooling  at  the  entrance  was  prevalent. 

In  computing  the  quality  of  the  hydrogen  at  any  axial  station  in 
the  tube,  it  was  assumed  that  an  unappreciable  amount  of  vapor  was  gen¬ 
erated  when  the  bulk  was  subcooled.  This  assumption  is  valid  for  low 
subcooling.  Thus,  quality  was  assumed  to  appear  in  the  tube  only  when 
saturation  conditions  were  reached.  A  saturation  length  I  has  been 
defined  as  the  axial  distance  from  the  entrance  of  the  tube  to  the  loca¬ 
tion  of  saturation  conditions.  In  computing  I,  the  bulk  ten^jerature 
and  pressure  profiles  are  required.  The  pressure  profiles  were  obtained 
from  the  data,  and  the  bulk  temperature  profiles  were  obtained  by 
computation: 


I  (BU 

Where  Cp  is  evaluated  at  (T^^l+i  + 

To  each  pressure  thei'e  corresponds  a  saturation  temperature;  thusj 
I  may  be  computed  by  locating  the  first  station  where  —  ^b,L,sat 

(TbjL,sat  ^S’Sed  on  p^) ,  retiring  one  station,  computing  Al,  and  add¬ 
ing  the  station  length  (see  sketch) . 


PL,sat-l  ^,sat 


sat-1 


t 

Length  =  Lg^t.i 

The  pressure  drop,  to  obtain  Al,  was  computed  from 

Ilj,sat  “  sat-1  "  ^ 

where 

_  ^  _  %,3at-l  ~  sat 
^  I'sat  “  ^sat-l 


(B2) 


(B2),.  and  a»  fuaotion.  ta- 

deseriBe  Pgat  ®-®  *  faction  of  rb,@a±>  an,  iterative  proeeaB  yields 

%8af  Then, 


.,  Tb^sat  ~  Tb,i,^aat-1 

(^) 


and 

^  =  Lsat-1  +  ^T-  (B3) 

Reference  14  presents  an  experimental  study  of  the  void  fraction 
(flow  area  occupied  hy  vapor)  for  water  at  2000  pounds  per  square  inch, 
which  Is  pertinent  to  the  definition  of  I .  The  attenuation  of  gamma 
rays  was  the  experimental  technique  used  to  measure  the  void  distrlhutlon. 
Appreciable  quantities  of  vapor  were  measured  for  subcooled  boiling. 

Heat  flux  and  the  model  of  the  two-phase  flow  (whether  stratified  or  dif¬ 
fuse)  had  a  marked  influence  on  correlating  the  void  fraction  with  the 
quality. 

While  these  water  data  cannot  be  applied  directly  to  boiling  hydro¬ 
gen,  they  do  point  out  the  limitations  of  the  assumption  used  in  computing 
hydrogen  quality.  An  estimated  comparison  of  the  actual  quality  distribu¬ 
tion  through  the  tube  with  the  computed  distribution  might  look  like  the 
following  sketch: 


For  the  data  reported  herein,  the  subcooling  was  moderate  and  the  vapor 
generated  to  I  was  small.  The  solid  curve  is  not  a  straight  line  be¬ 
cause  the  pressure  drop  through  the  tube  produced  a  superheat  condition, 
which  is  favorable  for  vapor  production. 

There  were  a  number  of  experimental  runs  in  which  saturation  condi¬ 
tions  occurred  before  the  fluid  reached  the  heated  portion  of  the  test 
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section.  For  these  cases,  a  negative  I  with  respect  to  the  test- 
section  inlet  was  considered.  These  cases  were  computed  separately  from 
those  in  which  I  was  positive  and  will  he  treated  in  the  discussion  of 
quality. 

It  is  speculated  that  the  analysis  used  is  relatively  accurate  in 
the  raidportlons  of  the  tube,  the  portion  used  for  the  heat-transfer  re¬ 
sults.  The  quality  in  the  suhcooled  entrance  region  and  in  the  exit 
section,  where  a  vapor  layer  and  end  losses  reduce  the  production  of  ad¬ 
ditional  vapor,  will  not  be  predicted  accurately. 


Quality 


The  quality  downstream  of  the  saturation  length  is  computed  from 


X 


(B4a) 


For  the  experlmentEil  runs  in  which  saturated  conditions  existed  up¬ 
stream  of  the  heated  test  section,  I,  in  (B4a),  was  taken  as  zero,  and 
a  mass  fraction  of  vapor  m  was  added  to  give  the  local  quality: 

+  m  ( B4b ) 


The  method  for  computing  m  can  be  best  explained  by  referring  to 
the  following  sketch: 


Approach  section 


Flow  IJ 

P. 

— 

- 

1 

Pi 

^  1  Pin^ 

Saturation 

^ - Heated  - m- 

portion 

point 


'Exit 


^sat’^sat  “  "^1 


where  T^^  is  the  bulk  temperatiu’e  measurement  at  the  entrance  to  the 
test  section,  and  p-|_  Is  the  corresponding  measured  pressure.  From 

the  pressure-drop  data,  a  station  in  the  unheated  tube  is  determined 
where  saturation  conditions  exist.  The  bulk  temperature  of  the  fluid  is 
assumed  constant  throughout  the  unheated  section  up  to  this  point. 
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OTie  enthalpy  change  frcaa  the  saturation  point  to  the  heated  tube 
Inlet  caused  by  evaporation  is  approximately  the  difference  in  the  satu¬ 
rated  liquid  enthalpies  between  these  stations.  Since  the  properties 
change  little  in  the  approach  section,  one  may  determine  the  quality  at 
the  entrance  by 


m  -  T.  (BS) 

A 

where  Cp  is  the  specific  heat  evaluated  at  (Tgat  Tin)/2. 

Tl  =  =  bulk  temperature,  and  T^jj  =  saturated  bulk  temperature  at 

inlet  (based  on  pressure).  There  were  no  runs  where  saturation  occurred 
before  the  approach  section. 

Density  and  Velocity 

The  density  at  any  axial  station  can  be  computed  from  knowledge  of 
the  quality? 


Pb  = 


1  -  X 


(B6) 


^g,sat  *^i,sat 


The  local  velocity  may  now  be  con^juted  from  one -dimensional 
continuity: 


^av 


w 

PpA 


(B7) 


Estimation  of  Heat  Losses 

Since  no  experimental  method  was  used  in  measuring  the  heat  losses, 
they  had  to  be  estimated.  To  facilitate  the  calculation,  only  very 
simplified  geometries  were  used.  It  was  found  that,  among  all  the  pos¬ 
sible  causes  for  heat  loss,  the  major  ones  were  radiation  loss,  heat 
loss  through  the  electrical  popper  flanges,  and  heat  loss  through  copper 
voltage  taps. 

Radiation  loss.  -  The  radiation  from  the  body  at  Tq  to  an  envelop¬ 
ing  body  at  Tcan  is  the  radiation  heat  loss; 
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where  the  temperature  of  the  stainless-steel  vacuum  container  is 

taken  as  500°  R,  the  radius  Tcan  inches^  and  the  emlsslvity 

®can  0.5. 

Conduction  loss  through  flanges.  -  A  first  approximation  to  the  • 
axial  conduction  through  one  end  of  a  hollow  cylinder  with  power  gen¬ 
erated  in  the  cylinder  is  as  follows: 


-  2 


"(^o  -  -  T4) 

L3  -  L4 
2 


(B9) 


where  Tq  is  0.1875  and  0..313  inch  and  r^^  is  0*1565  and  0.248  inch  for 
the  0.375-  and  the  0.625- inch  tubes,  respectively;  the  thermal  conduc¬ 
tivity  of  Inconel  or  stainless  steel  i.s  approximately  2X10"^ 

Btu/(in. )(sec)(°R) j  and  L3  -  L4  is  0.5  inch. 

Conduction  loss  through  voltage  taps.  -  Conduction  through  a  slim 
rod  is  obtained  as  follows : 


Qv  =  (bio) 

where  the  number  of  voltage  taps  ny  is  8,  the  radius  of  the  voltage- 
tap  wire  rjj  is  l/64  inch,  the  thermsil  conductivity  of  copper  /C(^^ 

is  4.63X10"'^  Btu/(in. ) (sec) (°R) ,  the  temperature  of  the  junction  strip 
for  instrumentation  within  the  vacuum  jacket  Tg  =  (Tq/2)  +  250°  R,  and 

the  distance  between  tube  and  instrument  junction  strip  Lj  -  Lg  is  2 
inches . 

In  general,  the  averaged  outside- wall  temperature  increased  with 
total  heat  input  Qtot>  term  Tg  -  T4  shows  the  same  trend. 

The  following  table  summarizes  the  estimates  of  various  forms  of  heat 
loss  from  the  Inconel  tube  as  a  function  of  total  heat  input  ^ot* 


^■Negative  sign  denotes  flow  of  heat  into  test  section. 

The  tabulated  values  of  heat  losses  indicate  that  data  corrections 
for  heat  loss  are  not  warranted. 
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PARA- ORTHO  CONVERSION 

Because  of  the  strange  shapes  of  vail  teiaperature  profiles  encoun¬ 
tered  In  the  hoiling  of  llc[uid  hydrogen,  q,uestions  were  raised  as  to 
the  possihle  effect  on  the  heat  transfer  of  the  conversion  of  para- 
hydrogen  (p-Hg)  to  ortho-hydrogen  (o-Hg).  This  speculation  was  "based 
upon  the  following  considerations:  ~ 

(1)  The  liquid  hydrogen  used  contained  about  95  percent  while 

the  equilibrium  concentrations  of  ^-Hg  are  50  percent  at  138°  R  and  25 
perOenL  aL  540°  H. 

(2)  The  heat  of  reaction  for  para- ortho  conversion  is  known  to  be 
comparable  to  the  latent  heat  of  vaporization  at  or  near  the  boiling 
point.  In  addition,  the  specific  heat  and  thermal  conductivities  of  the 
two  modifications  of  hydrogen  are  quite  different  when  the  temperature 
is  below  650°  R. 


Mechanism  of  Para-Ortho  Conversion 

In  order  to  estimate  the  effect  of  this  para-ortho  conversion,  it 
is  necessary  first  to  examine  the  associated  mechanisms  and  then  to 
determine  which  one  may  apply  to  the  experimental  system. 

The  para-ortho  conversion  may  take  place  in  the  following  forms : 

(1)  Homogeneous  reactions: 

(a)  Direct  transition:  A  very  slow  reaction,  about  one  tran¬ 
sition  every  300  years. 

(b)  Collision  of  molecules:  Half-life  time  is  3  years. 

(c)  Collision  with  hydrogen  atoms  (h) :  The  H  atoms  may  be 
generated  by  photons,  electric  charges,  alpha  rays,  and  gamma  rays. 
But  it  is  also  a  very  unfavorable  reaction  for  ortho  to  para  con¬ 
version  at  138°  R  (see  ref.  15).  Since  no  radiation  sovirce  exists 
in  the  hydrogen  heat-transfer  apparatus,  this  reaction  mechanism 
is  eliminated. 

(2)  Heterogeneous  reactions: 

(a)  Paramagnetic  mechanism:  The  Hg  molecule  is  absorbed  and 
the  spin  is  reversed  (ref.  16).  Catalysts  for  this  process  are 
paramagnetic  oxides,  charcoal  with  traces  of  in^iurities,  and  others. 
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(b)  Chemisorptions  Hg  molecttLss  are  oheraiserbed  ^on  metal  films 
and  activated  (ref*  16).  She  bonds  between  H  atoms  are  brb^&  and 
recombined.  Thus,  a  good  catalyst  must  have  a. strong  chemisorption 
property  toward  H2.  Among  the  active  chemisorption  materials  men¬ 
tioned  in  reference  17  are  tungsten,  nickel,  tantalum,;  and  tellurium. 
Thus,  it  follows  that  stainless  steel  and  Inconel  should  be  good 
catalysts.  Since  the  tube  wall  is  the  hottest  spot  in  the  appara¬ 
tus,  the  para- ortho  conversion  should  take  place  at  the  stainless- 
steel  tube  wall  if  there  is  any  such  reaction. 

Based  upon  the  preceding  survey,  it  appeexs  that  investigation  about  the 
possibility  of  para- ortho  conversion  should  be  directed  toward  the  reac¬ 
tion  at  the  wall. 


Conversion  Rate  in  Heterogeneous  Reactions 
The  mechanism  for  conversion  can  be  assumed  to  occur  as  follows: 

(1)  Diffusion  of  £-H2  through  a  boundary  layer  to  the  wall 

(2)  £-H2  to  0-H2  on  the  catalytic  wall 

(3)  Diffusion  of  0-H2  out  through  the  film 

By  assuming  the  para-ortlio  conversion  to  be  a  first-order  reaction. 


"^l^ortho,w  ^2^para,w 


^  dCpara 
-  »w  <ay 


(Cl) 


where  k^^  and  k2  are  absolute  reaction- rate  constants,  9  is  the  dif¬ 
fusion  coefficient,  and  y  is  the  distance  from  the  wall.  With  the  as¬ 
sumptions  of  a  lineeir  concentration  profile  across  the  boundary  layer  of 
thickness  &  and  it  can  be  shovm  that 


9  = 


^  Kkj^Cpaj.a^'b _ 

(1  +  K,^)ki  +  ^ 


(C2) 


where  9  is  reaction  rate  of  moles  per  unit  area  per  unit  time;  b  re¬ 
fers  to  bulk  temperature;  Keq  is  the  equilibrium  constant  at  wall  tem¬ 
perature,  Cortho,w/°para,w^  and  Kg^  =  3  for  T  >  350°  R. 


If  the  analogy  between  heat  transfer,  mass  transfer,  and  momentum 
transfer  is  used  (ref.  9),  it  can  further  be  shown  that 


.24 


^overall  “  thermal 


1  + 


(Kgq^V^c^)  Cpara 


vnij 


(l  +  KecL)ki  +  h^hermal  o^p 


or 


h  -  h  4  ^ 

overall  “  thermal  ^ 


(C3) 


where  h^yg^all  the  overall  heat- transfer  coefficient,  h+.hQTTpni  is 
thaL  for  the  case  without  para- ortho  conversion,  and  f^  =  ^p/|i  is  a 
constant  (f^  =  1.28  to  1.37) . 

However,  since  in  most  practical  cases  the  wall  temperature  is  near 
or  above  540°  R  where  M  is  very  small,  the  contribution  to  the  heat 
transfer  due  to  para- ortho  conversion  must  be  negligible.  !Phus,  even  if 
there  is  some  conversion,  the  only  effects  on  heat  transfer  will  be  re¬ 
lated  to  changes  in  the  magnitude  of  Cp  and  transport  properties.  As 

to  Tw  <  540°  R,  AH  is  not  zero,  and  equation  (C3)  may  still  give  some 
information  about  this  additional  heat  transfer  due  to  conversion. 


Estimation  of  Reaction  Rate 

To  weigh  the  effects  of  the  secondary  factors  such  as  differences 
in  thermal  conductivities  and  specific  heats,  it  is  instructive  to  esti¬ 
mate  the  reaction  rate.  Equation  (C2)  may  be  used  if  is  known. 

However,  because  of  the  lack  of  data,  k]_  is  not  readily  available.  The 
only  information  that  is  useful  is  an  experimental  determination  of  half- 
life  time  T  for  para  concentration  at  a  nickel  wire  with  surface  area 
of  1/2  square  centimeter  and  the  gas  volume  at  773°  K  and  p  =  100  milli¬ 
meters  of  mercury.  These  half-life  data  will  be  used  to  estimate  reac¬ 
tion  rates.  From  reference  9, 


Ten^. , 

T,  °K 

Half-life, 
T,  sec 

323 

5000 

398 

2000 

423 

800 

448 

500 

473 

150 

It  was  further  shown  in  reference  9  that 
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Tp  =  constant  (C4 

This  information  Is,  of  course,  far  from  acLeq.aate.  But  to  get  a  first- 
order  estimate,  this  is  better  than  nothing.  Further,  it  should  be 
noticed  that,  since  no  information  was  available  as  to  the  thickness  of 
film,  and  so  forth,  application  of  these  data  also  implies  the  assun^- 
tion  that  diffusion  rate  in  this  experiment  was  about  the  same  as  the 
one  in  the  present  system. 

In  reference  9, 

dCf 

•*  _  .v/r'.  .  n  ^  (nr 

dt  '  '  ^eq^ 

where  refers  to  concentration  at  time  t,  and  k  refers  to  reac¬ 
tion  rate.  In  terms  of  moles  converted  per  unit  time  and  area,  then, 

^  ~  ~  (^t,para  ‘  Ceq,para)  ^w 

=  (C6 


where 


_  kV  _  In  2  /V\ 
^  "  A  "  T  [a) 


But,  for  the  experimental  system, 


V  =  500  cc,  A  =  0.5  sq  cm 


Therefore, 


k'  «  - 


(in  2)  ^ 

0.5  693 


cm/ sec 


Suppose 

T^^  =  852°  R,  T  =  150  sec,,  p  =  100  mm 
At  conditions  where  p  -  800  millimeters,  using  equation  (C4),  T  is 


T  =  150 


/800\^ 


=  344  sec 


S6 


Further,  at 


then 


Ty  a  852^  R,  p  s  800  nni 


°t,para  =  ^  £S£2  ^^0“ 


n 

'^para^eq 


0.25  -E- 
RT 


Thus, 


«  -  693  ,  . 1  1^\ 

^  ~  344  ^  ^  ^760  82.05  852/ 

a  4.1X10"^  g-mole/(cm^) (sec)  =0.6  Ib/Csq,  ft) (hr) 

At  a  wall  temperature  of  500°  R,  ^  due  to  para-ortho  conversion 
will  he  small,  AP  large,  and  01  small j  therefore. 


^thermal  “  '^overall 

Knowing  01,  the  total  weight  converted  at  a  specific  flow  rate  can 
he  estimated.  The  total  conversion  rate  in  a  12-inch-long,  l/2- inch- 
ins  ide-diaraeter  tube  becomes 


^tot  =  =  0.0785  Ib/hr 


For  a  flow  rate  of  0.1  pound  per  second,  the  weight  percent  converted  is 


0.0785 
3600(0. l) 


w  0.0225^ 


Thus,  para- ortho  conversion  should  have  no  effect  on  the  experimental 
data. 


Experimental  Results 

Since  the  para- ortho  conversion  is  endothermic,  it  should  help  the 
heat  transfer;  and  thus  the  stainless-steel  tube  should  be  expected  to 
show  a  lower  wall  temperature  than  a  gold-plated  tube.  But  the  experi¬ 
mental  data  showed  the  gold-plated  tube  to  give  a  lower  wall  tonperature. 
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CS-4  back 


Therefore ,  the 
Is  more  likely 
conversion. 


deviation  between  the  two  prof lies  of  wall  ten^si^atnre 
to  be  due  to  experimental  error  than  to  para“6fth|0 


10 

r- 

o> 

I 


Conclusions 

(1)  The  para- ortho  conversion  should  take  place  on  catalytic  tube 
walls . 

(2)  The  reaction  rate  was  estimated  to  be  negligible  in  the  present 
system.  Further,  because  of  low  heat  of  reaction,  the  effect  on  heat 
transfer,  if  any,  should  be  only  secondary. 

(3)  Experimental  results  with  different  surface  materials  showed 
no  effect  of  the  presence  or  absence  of  a  catalytic  surface. 


f 


i 
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APPENDIX  D 


MAKPUffiLLI  PARAMETER  Xtt 
Introduction 

One  of  the  varlahles  used  in  correlating  the  data  is  the  Martinelli 
two-phase  parameter  Xtt*  two-phase  pipe-flow  model  was  composed 

of  a  cylindrical  gaseous  core  and  a  liq.uid  annulus  around  it  flowing 
along  the  wall  isothermally . '  The  general  form  of  the  two-phase  pressure- 
drop  equation  is 

The  parameter  T  must  he  determined  experimentally  and  is  a  cor¬ 
rection  coefficient  for  the  hydraulic  diameter  of  annular  flow  and  the 
velocity  of  one  phase  relative  to  the  other.  Martinelli  interprets  the 
Xtt  parameter  as  an  index  of  the  wall  and  two-phase  shear  resistance  of 

the  two-phase  flow. 

There  are  two  differences  in  the  two-phase  hydrogen  flow  model  as 
compared  with  Martinelli 's  original  model.  First,  heat  transfer  is 
present,  whereas  Martinelli  considered  only  an  Isothermal  case.  The 
presence  of  heat  transfer  would  affect  the  magnitude  but  not  the  defini¬ 
tion  of  the  friction  factor.  (Martinelli  used  the  Blasius  equation.) 

Second,  the  hydrogen  gas  phase  is  adjacent  to  the  wall,  and  the 
liquid  phase  makes  up  the  core.  This  is  an  inversion  in  the  position 
of  the  phases  in  the  pipe  as  specified  by  Martinelli.  However,  it  is 
shown  that  the  same  two-phase  parameter  Xtt  evolves  from  the  subse¬ 
quent  analysis,  which  is  quite  similar  to  that  of  Martinelli. 


Martinelli  Two-Phase  Analysis 
The  principal  assumptions  used  in  reference  1  are: 

(a)  Gas-phase  Ap/AL  equals  liquid-phase 

(B)  Volume  of  gas  plus  volume  of  liquid  equals  volume  of  pipe  at 
any  Instant. 

The  pressure  drops  (gas  and  liquid)  can  be  expressed  by  the  Weisbaoh 
equation  (ref.  l) : 


E-976 


S-976 


89 


^2 


-  „3 

Pi  ^ 

Dj  2gc 


Pg 

t  —2, 

Dg 


(D2) 


Note  that  the  Ap  of  the  two-phase  flow  is  greater  than  that  of  either 
single  phase  because,  in  addition  to  wall  friction,  there  exists  a  vary¬ 
ing  gas -liquid  Interface. 

Again  recall  that  Martinelli  considered  a  two-phase  system  in  which 
the  liquid  was  an  annulus  flow  adjacent  to  the  wall,  the  gas  phase  con¬ 
stituting  the  core.  Defining  the  hydraulic  diameter  for  cylindrical 
flow  as 


£  d2  =  A 
4 


where  A  Is  the  area  encompassed  by  the  fluid,  Martinelli  wrote  expres¬ 
sions  for  the  effective  hydraulic  diameter  of  each  phase' 


(D3) 


where  a  and  p  are  form  factors  and  in  effect  ratios  of  the  .actual 
cross-sectional  area  to  a  circular  cross  section.  For  his  flow  model, 
Martinelli  argued  that  P  «  1  and  a  »  1.  For  the  hydrogen  flow, 
a  1  and  p  »  1;  thus,  some  modification  to  his  presentation  is 
necessary. 

Assuming  a  =  1  and  P  1,  the  mean  fluid  velocities  may  be  com¬ 
puted  from 


4W7 


4wp 


(D4) 


so 


where  wj  eusid  Wg  are  flow  rates  of  the  .lig.uid  a»d  gas^  respeotively. 

Martlnelll  notes  that>  strictly  speajclng,  equations  (D4}  cannot  be  sub¬ 
stituted  In  equations  (IJ2).  However,  the  e3q)erlmentri  determination  of 
0  or  a  corrects  for  this  situation.  The  velocities  In  equations 
(D2)  should  Involve  the  velocity  of  one  phase  relative  to  the  other. 


The  Isothermal  friction  factors  and 
the  generalized  Blaslus  form: 


may  be  expressed  In 


(Ds; 


where  Kj  and  Kg  are  empirical  constants.  For  the  case  where  both 
phases  are  turbulent, 


Kg  =  K;  =  0.184 
n  »  m  =  0.2 


Similar  friction  factors  for  nonisothermal  conditions  may  be  defined, 
and  the  constants  change.  Using  equations  (D2)  and  substituting  equa¬ 
tions  (D3)  to  (D5), 


(D6a; 


(D6b; 
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Usiflg  assun©ttbn  (a) aijd  (D6b),  solve  for  the  ratio 
Dj/bg  and  substitute  the  values  for  m  and  ni 


Using  assun^tlon  (B)^ 

I  bfl  +  I  3b|L  =  I  b^ipeL  (b8) 


Df  +  pb2  «  D|ipg  (D9) 

Diameters  Dj  and  bg  are  unknowns.  Martihelli  solved  for  bg, 

which  was  the  diameter  of  the  cylindrical  gaseous  core.  For  the  hydro¬ 
gen  model j  Dj  Is  the  diameter  of  the  cylindrical  liquid  core.  Conse¬ 
quently,  an  expression  for  this  parameter  will  be  found  using  equation 
Cb7): 


This  is  similar  to  Martinelli ' s  answer  for  bg. 

Referring  to  equation  (D6a),  it  can  be  shown  that 


(DIO) 


(Dll) 


(D12) 
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9*hus , 


/PplpeV  ^ 

Wtp  ■  °l  ) 

■  ^).  W 


Using  equation  (Dll), 


f^)  =  (- 

Wtp 


/■  \-0.75/,,  \-0.083/„ \-0.416 

1  +  B0-25&1  ra  Pi 


Martlnelll  obtained  a  comparable  result: 


V^^p  Y^/g  W/  W/ 


+0.75/  v+0.083  -  ,+0.416 

1  /£s\ 


To  correlate  his  results,  Martinelli  sought  a  parameter  that  would, 
produce  a  for  various  fluids.  It  so  happened  that  the  multiplier  of 
a  in  the  previous  equation  did  the  Job.  For  simplicity,  this  was 
called  Xtt*  It  represents  the  shearing  force  in  the  two  phases. 

The  Xtt  parameter  was  found  more  convenient  to  compute  if  the  ex¬ 
ponents  were  changed  to  more  even  digits.  Raising  each  term  to  the  1.2 
power  does  not  change  the  "weight”  of  each  term;  thus,  X^t  written 

as  follows : 

(m  p) 

’^t  \nj  [^^1/  \pg/ 

/  \0.9/..  \0.1/oA0.5 


X  ]  fiis]  (n 

^  \^^l/  vs, 


In  this  report,  Xtt  used  rather  than  the  Inverse  form  as  would 
be  obtained  from  equation  (D14). 
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DESIGN  PROCEDUEE 

Ihere  are  many  design  methods  and  procedures  and  many  two-phase 
systems  in  which  the  data  of  this  report  may  be  applied;  however,  only 
one  system  will  be  considered  herein. 

In  applying  this  correlation  to  a  design  problem,  such  as  the 
channel  of  a  regeneratively  cooled  rocket  engine,  the  designer  should 
have  approximate  values  for  the  channel  geometry,  inlet  pressure  and 
temperature,  the  mass-flow  rate,  and  the  combustion-chamber  and  nozzle 
properties.  With  this  idea  of  channel  geometry,  the  designer  can  then 
calculate  the  heat-flux  distribution  with  or  wlthoi.it  radiation  by  a 
method  such  as  that  outlined  in  reference  18.  Methods  outlined  in  ap¬ 
pendix  B  will  be  employed,  without  further  reference,  to  obtain  some  of 
tlie  following  design  information. 


For  a  first  approximation,  assume  a  constant  pressure  profile,  p^^^. 

(a)  Determine  I  by  computing  the  bulk  temperature; 


^,L+1  =  '^,L  * 


%  ~  Ql\  M 
Lg  -  II/  -JAIp 


(El) 


where  (Qg  -  0,2) /(L2  -  is  the  heat  input  per  unit  length,  and  AL 
is  the  length  between  and  Lj^.  Since  Tb,5at  is  known  from  the 

pressure  (i.e., 

I  =  +  AZ  (B3) 


follows  immediately.  The  Z  estimated  in  this  manner  will  usually  be 
greater  than  the  actual  I,  but  is  a  good  approximation  for  low 
subcooling. 

(b)  Compute  x  using  X  evaluated  at  p  «  const,  a  Pin* 

(c)  The  average  density  at  each  station  is  then  estimated  using 
the  saturated  llq^uid  and  gaseous  densities.  The  average  density  prior 
to  saturation  ( Z  >  l)  is  obtained  from  the  bulk  temperatures,  previously 
computed  while  determining  Z. 

(d)  The  static  pressure  at  each  station  is  computed  using  one¬ 
dimensional  continuity  and  momentum  pressure  drop  (see  also  fig.  4): 
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w  a  puA  (02) 

-  ^  e«  _!L  /5U  (2) 

Age 

The  designer  must  now  Iterate  the  processes  of  finding  If 

Pl,  while  using  Pj;,  to  determine  Tb,L,sat  (®®®  ^^8*  8(a)),  until  the 
desired  accuracy  is  attained.  After  obtaining  this  information,  com¬ 
pute  Eq.uation  (6)  will  then  yield  the  ratio  .Nu5alc,f/NUexp,f • 

Compute  WuQaic,f  i-ising 

Nucalc,f  =  0*023  Re°-®Pr?‘^  (3) 

and  local  properties  evaluated  at  film  temperature  (arithmetic  mean  be¬ 
tween  bulk  and  v/all).  Transport  properties  are  available  in  literature 
from  the  National  Bureau  of  Standards  (refs.  5  to  7).  The  local  heat- 
transfer  coefficient  htp,L  now  becomes  available  through  Nuexp,f  (eq,. 
(6)).  Compute  the  coolant  wall  temperature  accounting  for  the  material 
and  find  the  heat-flux  distribution  into  the  fluid.  Since  input  heat  and 
absorbed  heat  will  generally  not  agree,  an  iteration  of  the  entire  bulk 
cooling  process  will  be  required  to  balance  the  gas-side  and  liquid- 
side  heat  fluxes.  One  must  keep  in  mind  that  the  htp,L  computed  will 
probably  be  conservative,  since  a  rocket  channel  heats  only  a  portion 
of  the  wetted  perimeter,  while  in  this  experiment  the  entire  wetted 
perimeter  was  heated. 

To  recapitulate,  using  figures  and  formulas  (see  text  and  appendix 
B  for  source  of  equations) : 

(1)  The  designer  must  know  something  of  the  geometry  and  condi¬ 
tions  to  be  met.  Consider  a  single- pass  channel.* 


Ap,  is  not  necessarily  >  A^+p 
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(8).  The  heat-flux  distrihution  must  be  .estimated  (see  ref . '18) : 


(3)  The  saturation  length  may  be  estimated  from 


.  L  * 

1  ^ 

(El) 

^L2  -  Lij 

wGp 

^sat-l  + 

(B3) 

For  constant  heat  Input  per  unit  length  or  saturation  prior  to  heating; 
see  appendix  B. 

(4)  The  local  q,uality  may  now  be  computed: 


(E3) 


and,  for  constant  heat  input  per  unit  length, 


X  =  ^ 

LfcW  \  ^  / 


LfcW 

(5)  The  local  two-phase  density  is  given  by 

1 


X  +  1  -  X 
Pg, sat  Pi ,  sat 


(B4a) 


(B6) 


(6)  The  local  velocity  may  now, be  estimated  from  the  one- dimensional 
continuity  equation: 


w 

PfaA 


(B7) 


(7)  The  local  pressure  drop  Is  estimated  using  equation  (2); 


-Ap=_iLdja  (2) 

ASc 

(8)  Iterate  the  preceding  steps  to  obtain  better  approximations 
for  the  flow  parameters. 

(9)  The  approximated  pressure  drops  niay  appear  as  follows: 


in 


Computed  from 
iteration  of  entire 
process  (not 
necessarily  lower) 


Assumed  for  first 
approximation 


Computed  from  first 
approximation 


Con5)uted  from  second 
approximation 


^  L 


(lO)  Calculate  from 


1 

^t,f 


(D16) 


(11)  Estimate  the  Nusselt  number  ratio  by  using  equation  (6): 


_ ^^calc>f _ 

"'^expjf  =  0.611  +  1.93  Xtt.f 


(6) 
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(12)  Compute  the  loceJ.  Nusselt  numher  using  the  modified-  Dittus- 
Boelter  equations 


Nucalc,f  =  0*023  Re°'®Pr?*^ 

vhere 


Re  =f£i2li£lP 


and 


(3) 

(4) 


_ 1 _ 

x/pf  +  (1  -  x)/p2 


(5) 


(13)  Using  steps  (12)  and  (ll),  compute  the  local  two- phase  heat- 
transfer  coefficient: 


1  F^^exp,f^ 


^tp,L  -pi  “ 


f  /L 


(E4) 


(14)  Compute  the  liquid-side  wall  temperature  Tw,  L  accounting 
for  the  wall  thickness : 

Other  forms  of  computing  heat  transfer  through  the  wall  may  "be  more 
desirable  and  should  then  be  used. 

(is)  The  heat  absorbed  by  the  liquid  may  now  be  estimated r 

^•Jabsorbed  by  liquid^L  “  Hp,L^Ar,  L  ”  ^b,L) 

(16)  Compare  <labsorbed,L  ^7  liquid  with  the  desired  qp, 
and  adjust  parameters  used  accordingly  and  iterate. 
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SUmARY  01'  HYOAOaeN  U£AT-'I1UN8PSR  DAHIA  PQR  0.37&*XtIQH«0Um0BrPXWT£A- iNCODifilj  TUBS 
[q  •  flow;  ^  •  hoat  fiux)  ft  «  ni?.0s-fiow  j'aWt] 


Axial  distance 
from  tube  Inlet, 
L, . 

In. 


Local  static 

Local 

Average 

Local 

Local  bulk 

mslde  wall 

pressure. 

velocity. 

density. 

quality, 

temperature. 

temperature. 

Ib/sq  In.  aba 

rt/aee  ■ 

“av' 

Ib/ou  ft 

X, 

$  vapor 
by  mass 

Tb* 

OR 

^W,li 

OR 

Heat-transfer 

coefficient, 

A, 

Btu/sq  in.-8eo-®R 


0.177  Ib/scc)  Q  «  2.77  Btu/aecj  q  «  0,235  Bbu/(8q  ln.)(3ec) 


0 

44.5 

249.0 

.004 

44.5 

192.1 

.010 

44,4 

229.6 

.021 

44.3 

255,9 

.032 

44.2 

2f>6.7 

.043 

44.0 

277,6 

.055 

43.9 

256.1 

.0G2 

43.7 

252.7 

.069 

43,5 

219.9 

.077 

43.3 

2G0.2 

.064 

43.0 

166,1 

.080 

42.9 

261.0 

Run  lH-3:  ft  «  0.136  Ib/ooc;  Q  «  2.71  Btu/aec;  q  «■  0.230  Btu/(8q  lr.)(oec) 


249.3 

0.00112 

236.0 

.00116 

.70  Btu/occj  q  ■  0,22'J  lilu/(oq  lr.)(occ) 


0.0*;.'.  in/jw<*j  Q  »  2.73  Biu/aec}  q  -  0.231  Utu/(aq  ln.)(oce) 


0.305  Blii/(sq 


05.6 

- 1 

2.G3 

0.026 

42.1 

203,0 

o.ooloe 

73.9 

•^.35 

.035 

42.0 

240. C 

.00115 

‘»3.1 

2.09 

.041. 

42.0 

263,0 

.00104 

102. y 

i.nu 

.06!, 

4 1 ,  Cl 

270.4 

.00100 

123.1 

1.41 

.065 

41.6 

260,7 

.00102 

144.fi 

1.20 

.105 

41,4 

264.2 

.00103 

166.0 

1.03 

•  l25 

41,2 

253.2 

.00119 

lv4.9 

.93/ 

.130 

41.0 

215.6 

.00131 

202.3 

.yi'V 

.151 

40,6 

212.2 

.00134 

219.0 

.  iC4 

40.6 

102. D 

.00151 

240.  i 

.  /22 

.177 

40,4 

07.2 

.00402 

.  (9i3 

.lu5 

40.2 

2C0.& 

.00104 
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TAELS  I.  -  CONTINUSD.  NUMMARY  OF  HYI^OQKN  HBAT-TRANSm  DATA  FOR  O.STS-IHW-OtrraiDE.DIAMGm  XNpp^  TOBB 
[|q  «  heat  flow)  q  m  heat  fluxj  W  •  masB-flow  rate. 3 


Axial  distance 
from  tute  Inlet, 
h, 

In. 


Looal  static 
pressure, 

Ib/aq  In.  abs 


Local  Average  Looal 
velocity,  density,  quality, 
u,  X, 

ft/oBc  ft  i,  vapor 

by  mass 


Looal  bulk  Inside  wall 
temperature,  temperature, 


Meat**  transfer 
coefficient, 
h, 

Btu/sq  ln,-aec-°R 


0.125  Ib/scc;  Q  »  4.61  Btu/sec;  q  m  0.590  Btu/(5q  ln.)(sec) 


Run  lO^b:  ii  M  0.088  Ib/ocoj  Q  ■■  4.C3  fitu/so 


0.392  Btu/(cq  ln.)(Dec) 


59.4 

2.78 

0.025 

43.1 

243.1 

71.6 

2.31 

.041 

43.0 

357.6 

65.4 

1.94 

.058 

43.0 

370.3 

115.3 

1.43 

.095 

42.9 

358.  B 

146.1 

1.15 

.132 

42.7 

3S1.0 

178.0 

.  925 

.168 

42.5 

346.5 

214.8 

.770 

.205 

42.2 

313.7 

240.9 

.686 

1  .230 

42.0 

30C.2 

269.4 

.614 

.252 

41.7 

291.9 

299.  / 

.552 

.2/5 

4J  .  4 

2/0.4 

55G.4 

.492 

.297 

40.9 

190.4 

361.2 

.450 

j  .311 

40.6 

242.6 

O.r/0  Ih/nccj  Q  m  7.0ft  Dlu/occ;  q  •  0.597  Blu/(oq  in.)(DCc) 


0.055 

!  57.9 

.64 

!  57.6 

1.3C 

.57.2 

2.ar. 

56.1 

4.3C 

54.9 

5.86 

!  53.5 

7.40 

1  51.7 

8.44 

50.2 

9.40 

48.1 

10.40 

46. G 

11.40 

44.2 

12.00 

42.4 

Run  20-2:  w  »  0.15C  Ib/oocj  Q  «  7.12  Blu/soc’j  q  •»  0.C05  Dlu/(8q  ln.)(8ec) 


Q  «*  7. 15  Btu/ 


oq  ln.)( 


I  I  \  -.b08Tffll®8.  aiMMAfty  b?  KYBRObKH  iaAT-*!Ujf!3!®  DATA'PPR  0.37Sini(9|-6U®iDE-DtAiffiTBR  IMM 


[q  «  heat  flow)  (j  ■  heat  flux)  A  ■  maaa-fjow  rate,] 


Axial  distance 
from  tube  inlet, 

J' 

in. 

Local  static 
proBsure, 

Ib/aq^in.  ahs 

Local 

velocity, 

u, 

ft/sao  . 

Average 

denelty, 

Pav' 

Ib/ou  ft 

Local 

quality, 

4  vapor 
by  mass 

Local  bulk 
temperature, 

Inside  wall  . 
temperature, 

°R 

Heat-transfer 

ooefficlent, 

h, 

Btu/sq  ln,-8eo-®R 

Run  20-4 i  w 

•  0.079  Ib/aeo)  Q  -  7. 

17  Btu/sec|  q  ••  0.608  Btu/(eq  ln.)(8eo 

■ 

0.055 

37.4 

42.3 

3.4B 

0.009 

43.0 

279.2 

0.00258 

.64 

37,0 

61.4 

2.40 

.037 

42.9 

545.4 

.00121 

1.36 

36.5 

63.3 

1.77 

.067 

42.8 

545.9 

.00121 

2.86 

35.4 

130.9 

1.13 

.131 

42.6 

513.0 

.00129 

4.36 

34.0 

161.9 

.809 

.194 

42.3 

502.7 

.00132 

5.86 

32.5 

237.1 

.621 

.257 

41.9 

502.7 

.00132 

7.40 

30.7 

300.3 

.490 

.320 

41.5 

579.1 

.00139 

6.44 

29.4 

347.5 

.424 

.362 

41.2 

462.2 

.00138 

6.40 

27.6 

400.8 

.367 

.399 

40.7 

453.5 

.00147 

10.40 

26.2  ' 

4S6.8 

.322 

.439 

40.3 

453.6 

.00147 

11.40 

24.1 

530.2 

.270 

.477 

39.7 

406.4 

.00166 

13.00 

22.  i 

OU£.0 

_ -^pp 

.499 

09. 

.00516 

Run  22-4 :  m 

>  0.151  Ib/aco)  Q  .  8. 

23  Btu/sec;  q  ■  0.697  Btu/'(8q  in.}(8ec 

0.05S 

65.0 

71.4 

3.9? 

0 

45.4 

259.4 

0.00326 

.64 

64.7 

72.3 

3.92 

0 

46.2 

530.6 

.00144 

1.36 

64.4 

73.4 

3.06 

0 

47.2 

565.9 

.00134 

2.86 

63. 4 

100.4 

2.82 

.034 

47,5 

533.7 

.00130 

4.36 

62.0 

134.5 

2.11 

.076 

47.3 

563.1 

.00135 

S.86 

60.2 

170.7 

1.66 

.J17 

47.0 

529.6 

.00144 

7.40 

S6.0 

210.9 

1.34 

.159 

46.7 

494.6 

.00156 

8.44 

56.3 

239.4 

1.18 

.186 

46.4 

480.. 1 

.00161 

8.40 

54.0 

?6fl.5 

1.06 

.2)1 

46.0 

462.2 

.00168 

10.40 

51.4 

301.9 

.930 

.236 

45.6 

444.6 

.00175 

11.40 

48.1 

341.3 

.030 

.259 

45.1 

410.6 

.00191 

12.00 

45.0 

368.9 

.766 

.273 

44.7 

260.2 

.00324 

Run  20-5:  w 

.  0.178  Ib/tccj  «  -  9. 

77  Btu/occj  q  »  0.828  Btu/(8q  ln.)(aec 

0.0S5 

66.4 

80.4 

4.14 

0 

42.3 

236.0 

0.00428 

.64 

66.1 

81.4 

4.09 

0 

43.5 

594.2 

.00150 

1.36 

65.6 

02.7 

4.02 

0 

44.4- 

650.0 

.00137 

2.06 

64.6 

05.4 

3.90 

0 

46.6 

639.2 

.00140 

4.36 

63.4 

106.7 

3.12 

.022 

47,5 

610.5 

.00145 

5.06 

61.7 

147.2 

2.26 

.064 

47.3 

601.4 

.00150 

7.40 

59.5 

191.0 

1.73 

.107 

46.9 

551,7 

.00164 

0.44 

57.7 

224.0 

1.49 

.135 

46. 6 

527.3 

.00172 

0.40 

55.0 

257.1 

1.29 

.160 

46.2 

479.5 

.00191 

10.40 

52.0 

’92.0 

1.14 

.186 

45.8 

456.6 

.00202 

11.40 

43.4 

135.4 

.992 

.211 

45.3 

400.2 

.00228 

12.00 

46.7 

356.0 

.900 

.225 

44.8 

188.5 

.00576 

Run  20-Ct  w 

•  O.j.32  Ib/occ;  Q  •  0. 

93  Btu/cccj  q  -  0.841  13tu/(nq  ln.)(occ) 

0.005 

57.9 

0.5 

4.10 

0 

43.0 

250.2 

0.00406 

.64 

57.6 

6... 5 

4.03 

0 

44.3 

631.0 

.00143 

1.36 

56.3 

62.0 

3.95 

0 

45.8 

653.2 

.00138 

2.86 

55.7 

97.x 

2.55 

.045 

46.5 

605.1 

.00150 

4.36 

54.3 

143.2 

1.73 

.099 

46.1 

566,3 

.00102 

5.06 

52.6 

131.7 

1.29 

.]a''>4 

45.8 

540.1 

.00168 

7.40 

50.1 

246.1 

1.01 

.210 

45.1 

500.6 

.00182 

8.44 

40.2 

266.2 

.666 

.246 

45.1 

507.6 

.00162 

3.40 

45.4 

330.2 

.750 

.279 

44.6 

460.0 

.00193 

10.40 

43.3 

376.6 

.656 

.313 

44.2 

477.4 

.00194 

11.40 

33.9 

436.3 

.568 

.344 

43.5 

432.6 

.00218 

12.00 

37.4 

460.6 

.516 

.363 

43.0 

196.0 

.00543 

Run  20-7:  w 

»  0.093  Ib/oec)  0-9. 

93  Btu/occj  q  >■  0.841  Bta/(oq  .In.) (see 

■ 

0.055 

49.3 

42.9 

4.05 

0 

44.0 

267.6 

0.00376 

.64 

49.0 

49.1 

3.53 

.000 

45.2 

653.2 

.00138 

1.36 

48.5 

73.0 

2.38 

.046 

45.1 

653.7 

.00138 

2.66 

47.4 

124.2 

1.40 

.125 

45.0 

577.7 

.00158 

4. 36 

4C.1 

177.7 

.975 

.202 

44.7 

541.9 

.00169 

5.66 

44.3 

235.4 

.737 

.278 

44.4 

559.5 

.00163 

7.40 

41.0 

301.6 

.S75 

.355 

43.9 

546.0 

.00168 

6.44 

33.3 

3S1.9 

.493 

.405 

43.5 

555.4 

.00164 

9.40 

36.9 

411.6 

.422 

.449 

42.9 

525.3 

.00174 

10.40 

34.6 

473.0 

.367 

.496 

42.5 

526.8 

.00174 

11.40 

31.1 

5C2.2 

.309 

.539 

41.6 

477.4 

.00193 

12.00 

26.1 

636.9 

.272 

.563 

40.8 

211.5 

.00493 

!i 
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Figure  1.  -  Schematic  diagram  of  liquid-hydrogen  heat-treuisfer  apparatus. 


SdieiBstlc  or  e^tpantiis  aad  Instrumeatatlan*  (All  dlmenBlooB 


Length;  L;  ln> 

(a)  0.375-Inoh-out8ide-dlaineter  Ineonel  tube. 
Flgttre  3.  -  'laical  Inalde-vall  temperatiire  dlstrlhutlons 


Length,  L,  in. 


(h)  0.625-Inch-outBide-diaiDeter  stainless- steel  tube. 
Figure  illustrates  vapor  hindii-jg. 

Figure  3.  -  Concluded.  Oypical  inside-wall  temperature 
d'Lstrihutions. 


ilet 


Figure  4.  -  Compa 
along  length  of 
Inconel  tube. 
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2  4  6  S  10 

Lengthy  Zj,  In- 

(b)  0- 625-Inch-out8lde-aiaineter  stainless-ateel  tube. 


Figure  5.  -  Concluded.  lyploal  local  values  of  heat-transfer 
coefficient  along  length  of  test  section. 


Nvisaelt  number  ratio. 


120 
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Figure  8.  -  ThermoaynBmlc  and  transport  data  for  liquid  and  gsiseciUB  hydrogen. 


Figure  8.  -  Contlnueii.  Thermodynamic  and  transport  data  for  liquid  and  gaseous  hydrogen. 
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(2)  Liquid  para-hydrogen  (refs.  5  and  6). 
(t))  Concluded.  Specific  heat. 


Figure  8.  -  Continued.  Thermodynamic  and 
transport  data  for  liquid  and  gaseous 
hydrogen. 


Figure  8.  -  Continued.  Thermodynamic  and  transport  data  for  liquid  and  gaseous 
hydrogen. 


(2)  GsiseouB  hydrogen  (refs.  5  and  6). 

(c)  Concluded.  Yiscosity. 

Figure  8.  -  Continued.  Theme dynctmlc  and  transport  data  for  liquid  and  gaseous  hydrogen- 
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(e)  Density- 

Figure  8.  -  ContlnueiZ.  IliejnBoaynaaii.o  euid  tronsport  data  for  liquid  and 
gaseous  hydrogen- 
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